Probing new physics in 5 ^ 07r decay 
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Abstract 

We analyze the rare decay mode B (pir, which is a pure penguin induced process, 
receiving dominant contribution from the electroweak penguins. Thus the standard 
I model branching ratio is expected to be very small, which makes it as a sensitive probe 

of new physics. Using QCD factorization approach, we find the branching ratio in 
I the standard model as Bv{B^ 4''^^) — 5 x 10^^. Exploring some of the beyond 

'nJ" \ standard model scenarios the branching ratio is found to be ~ O(10~^) in the minimal 

\ supersymmetric standard model (MSSM) with mass insertion approximation and in the 

\ extended technicolor model (TC2). The existence of an extra vector-like down quark 

(VLDQ) model predicts it to be ~ ©(lO"'^). The recent BaBar result on 4>tt^ 



p. 



^ , might be a strong indication of new physics effect present in the penguin induced 

D , process B ^ (pn. 



X 



PACS : 13.25.Hw, 12.60.Jv, 12.60.Nz 



1 Introduction 



The intensive search for physics beyond the standard model (SM) is performed now a days 
in various areas of particle physics. In this respect, the B meson system can also be used as a 
complementary probe to the search for new physics. The main objectives of the ongoing and 
the future B factory experiments are to explore in detail the origin of CP violation, to test 
the standard model at an unexpected level of precision and to look for possible existence of 
new physics effects. In the B experiments new physics beyond the SM may manifest itself in 
various ways, e.g., (i) decays which are expected to be rare in the SM and are found to have 
large branching ratios (ii) CP violating rate asymmetries which are expected to vanish or to 
be very small in the SM arc found to be significantly large (in) the discrepancy between the 
mixing induced CP asymmetry between various B decay processes which are dominated by 
a single decay amplitude with same weak phase, i.e., S^Ks ^-nd S^Ks ^tc. 

Thus the rare B meson decays are suggested to give good opportunities for discovering 
new physics beyond the SM. The discrepancy between the recently measured S^Ks ^-iid S^Ks 
[1, 2] has already given an indication of the possible existence of new physics in the B decay 
amplitudes (i.e., in the penguin induced B — > (t)Ks decay). 

In this paper, wc would like to explore the presence of new physics in another pure penguin 
induced decay mode B 07r. It proceeds through the quark level transition h dss, which 
is a flavor changing neutral current (FCNC) process at the one loop level. The interesting 
feature of this process is that it is dominated by electroweak penguin contribution. The QCD 
penguins should play a minor role for this transition since s and s quarks emerging from the 
gluons of the usual QCD penguin diagram form a color octet state and consequently cannot 
build up the (f) meson which is a ss color singlet state. Therefore the SM prediction for the 
branching ratio for this process is expected to be quite small. Thus as this mode is highly 
suppressed in the SM, it may serve as a good hunting ground to look for new physics beyond 
the SM. 

Recently, these decay modes have been searched for by the BaBar collaboration [3, 4]: 

Bt{B^ (j)n^) < 0.41 X 10"^ , 
Br(5° ^ 07r°) = (0.2+°;^ ± 0.1) x 10"^ , 

< (1.2±0.8) X 10-^ (1) 

On the theoretical side the decay mode B~ 07r^ has been studied recently in Ref. [5]. 
Using QCD factorization approach the branching ratio has been obtained in the SM and in 
the constrained minimal supersymmetric standard model. It has also been studied in Ref. 
[6] in the SM and R-parity violating sypersymmetric model and in the Refs. [7, 8] using the 
standard model approach. 

In this paper, we will first reanalyze this process again in the SM for the sake of complete- 
ness, using QCD factorization. It should also be noted that there is not good agreement of 
the SM predictions in the literature. We then consider the minimal supersymmetric model 
with mass insertion approximation. Basically we will confine ourselves to the case where the 
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new contributions to the 5 — > 07r will arise from the gluino mediated b — > dss process in- 
duced by the flavor mixing in the down s-quark sector. Next we will calculate the branching 
ratio in the framework of topcolor assisted technicolor model (TC2) and in the model with 
an extra vector like down quark (VLDQ). 



2 Standard Model contribution 



In the SM, the decay process B ^ (1)71 receives contribution from the quark level transition 
h dss, which is induced by the pure penguin diagram with dominant contributions coming 
from electroweak penguins. The effective Hamiltonian describing the decay b — > dss [9] is 
given as 

r [ 1 

^ W.E^A , (2) 

'- 1=3 



njSM 
^eff 



V2 



where q — u, c. O3, ■ ■ ■ ,0q and O7, ■ ■ ■ , Oiq are the standard model QCD and electroweak 
penguin operators, respectively. The values of the Wilson coefficients at the scale ~ m;, in 
the NDR scheme are given in Ref. [10] as 

C3 = 0.014 , C4 = -0.035 , C5 = 0.009 , Cq = -0.041 , 

C7 = -0.002Q; , Cg = 0.054Q; , Cg = -1.292q; , do = 0.263q; . (3) 

We use QCD factorization [9] to evaluate the hadronic matrix elements. In this method, the 
decay amplitude can be represented in the form 



{(f)n-\Oi\B-) = {(f>n-\Oi\B- 



/fact 



l + Er„< + C(AQCD/mb) 



(4) 



where (07r~|Oi|-B~)fact denotes the naive factorization result and Aqcd ~ 225 MeV, the 
strong interaction scale. The second and third terms in the square bracket represent higher 
order as and Aqcd/^t.^ corrections to hadronic matrix elements. 

In the heavy quark limit the decay amplitude for the B~ — > 07r~ process, arising from 
the penguin diagrams, is given as 



\SM, 



B- 



G 



V q=u,c 



«3 + «5 - 2 + Og) 



(5) 



where X is the factorized matrix element. The amplitude for B^ 07r° is related to 

B~ — >• 07r~ by 



(6) 
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Using the form factors and decay constants, defined as [12] 



{7:-{p^)\d-i^h\B-{pB)) 



— m: 



2 n 



(0(g,e)|s7^s|O) = /^m^e^ (7) 
we obtain the factorized matrix element X as 

X = (7r-(p.)|J7M(l-75)&|5-W)(0(g, 6)1^7^1 -75)s|0) 

= 2Ff-^ml) Um^ie-ps). (8) 

The coefficients a1,s which contain next to leading order (NLO) and hard scattering correc- 
tions are given as [7, 11] 



«3 








< 
















< 





(9) 



where = 3, is the number of colors and Cp — {N"^ — l)/2N. The parameters present in 
(9) are given as 



F4> 
ell 



-121n-^-18 + 4 + /f , 



/ dx g{x)(f)^{x) , 
Jo 

l-2x^ 
6- m x — Sin , 



1 — X 



B 



1 dz 



N Ff-^^(0)m| Jo z 



1 dx 



X 



dy 
y 



My) ■ 



(10) 



The light cone distribution amplitudes (LCDA's) at twist two order are given as 



2 2 

m^x 



(j>B{x) — Nbx^{1 — a;)^exp^ 

(t>%A^) = - ^) ' (11) 
where Nb is the normalization factor satisfying dx(l)B{x) = 1 and lob — 0.4 GeV. The 
branching ratio can be obtained using the formula 



!)7r-) 



Tb- 



2 



be 



STTirij 



\A{B- ^<P7T-)/{e-pB) 



(12) 



where k — Tgo / tb- and Pcm is the momentum of the outgoing particles in the B meson rest 
frame. 

For numerical evaluation we have used the following input parameters. The value of 
the form factor at zero recoil is taken as F^~*'^{0) = 0.33, and its value at = m| can 
be obtained using simple pole dominance ansatz [12] as F^^^iw?^) = 0.34. The values of 
the decay constants are as = 0.233 GeV, fs = 0.19 GeV, /7r=0.131 GeV, the particle 
masses and the lifetime of B mesons tb- — 1.674 ps, Tbo = 1.542 ps are taken from [13]. For 
the CKM matrix elements, we have used the Wolfcnstcin parameterization and have taken 
the values of the parameters A = 0.819 ± 0.040, A = 0.2237 ± 0.0033, p = 0.224 ± 0.039 
and 7] = 0.324 ± 0.039. With these input parameters, we obtain the branching ratio for 
B- 07r~ in the SM as 

Bt{B- (f)Tr-)\sM = (5.5±0.9) X 10"\ 

Br(5° 07r°)|5M = (2.5±0.4) X 10-^ (13) 

These predicted values are quite below the present experimental upper limit and the central 
value of 5° — > 07r° (1). It should be noted that our prediction is in agreement with [5], the 
shght difference is due to the difference in the used CKM parameters and formfactor. 



3 Supersymmetric contribution 

Now we study the decay process B — > (pir, in the minimal supersymmetric standard (MSSM) 
model with gluino contributions, because the chargino and charged Higgs contributions are 
expected to be suppressed by the small electroweak gauge couplings. Thus, the one loop 
contributions to the above mentioned decay process can be induced by s-quark and gluino 
penguin and box diagrams. These gluino mediated FCNC contributions are of the order 
of strong interaction strength, which may exceed the existing limits. It is customary to 
rotate the effects, so that they occur in s-quark propagator rather than in couplings and to 
parametrize them in terms of dimensionless parameters. Here we work in the usual mass 
insertion approximation [14, 15] where the flavor mixing i j in the down-type squarks 
associated with qs and q^ arc parametrized by with A, B = L, R and i,j as the 

generation indices. More explicitly (5^^)^ = (V^^MjV^)jj/m?, where Mj is the squared 
down squark mass matrix and rriq is the average squark mass. Vd is the matrix which 
diagonalizes the down quark mass matrix. 

The new effective AB — 1 Hamiltonian relevant for the B ^ (frn process arising from 
new penguin/box diagrams with gluino-squark in the loops is given as 

^susY ^ -^v,,V;, E [Cro. + CrO.] , (14) 

where Oi are the QCD penguin operators and the Cj^^ are the new Wilson coefficients. The 
operators are obtained from Oj by exchanging L ^ R. Thus including the new physics 
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contribution, one can write the total amplitude for S — > ^tt process as 

j^SUSY 



\SM 



1 + 



JO 



(15) 



where A'^^'^"^ is the new physics amplitude arising from minimal supersymmetric model 
and 9 is the relative phase between the SM and the new physics decay amplitudes. The 
corresponding branching ratio is given as 



Br = Br 



SM 



l + 2cos^|A 



SUSY I ASM 



SUSY I aSM\2 



(16) 



where Br*^^ is the SM branching ratio. To evaluate the amplitude in the MSSM, we have to 
first evaluate the Wilson coefficients at the h quark mass scale. At the leading order in mass 
insertion approximation the new Wilson coefficients corresponding to each of the operator 
at the scale ~ m ~ M^^^ are given as [15] 



cr 



AGFVt,V,*,ml 



{Si.) 



13 



13 



13 



13 



-Ib,{x) + ^B,{x)-^P,{x) + ^P,{x) 
^B,{x) + ^B2{x) - ^Pi{x) - ^P2{x) 
~B,{x) + ^-B^ix) + ^Pi(x) + ^P2{x) 



(17) 



The corresponding Ci are obtained from Cf^ by interchanging L ^ R. The functions appear 
in these expressions can be found in Ref. [15] and x = m^/m^. Because the parity of the 
vector meson (f) is opposite to that of B and tt mesons, which are pseudoscalars, the gluino 
loop effects appear as (^f^)^^ + {Srr)_^^- The Wilson coefficients at low energy C^^(iji), 

jjL ~ 0{mb) can be obtained from C^^{Mw) by using the Renormalization Group (RG) 
equation as discussed in Ref. [10] , as 



C(/x) = U5(/i,MH.)C(M, 



w 



where C is the 6x1 column vector of the Wilson coefficients and U^in, Mw) is the five-flavor 
6x6 evolution matrix. In the next-to-leading order (NLO), U5(/x, Mw) is given by 



U5(/i,MH.) 



(19) 



where IJ f\iJ,,Mw) is the leading order (LO) evolution matrix and J denotes the NLO 
corrections to the evolution. The explicit forms of U5(/i, Mw) and J are given in Ref. [10]. 
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For the numerical analysis, we fix the SUSY parameter as rriq — rrig — 500 GeV, 
as{Mw) = 0.118, asimj, = 4.4 GeV) = 0.221. Thus, the values of the Wilson coefficients 
evaluated at the b quark mass scale are given as (where the common factor ((^ll) has been 
factored out) 

C^^ = 0.025 , Cf ^ = -0.021 , C^^ = -0.003 , C^^ = -0.087 . (20) 

Now evaluating the matrix elements of the operators Os-q as done in Eqs. (5)- (9) for the SM 
i.e., replacing the SM Wilson coefficients 6*3-6 by their corresponding SUSY counterparts, 
we obtain the total amplitude in the MSSM as 

A(B- ^ (/.TT-) = + 2.525 ((^^3)^^ + {S'Is)rr)) ■ (21) 

The constraint on (5^3 )ll is obtained from — mixing, however for {Sf^)Rji, it is not 
very stringent. We use the conservative hmit for mg=500 GeV and x = 1 as (^^3)^^^ < 



\^c{6i^)1^\ ~ 9-8 X 10 ^ [15], {Sf^)LL ~ {Si3)rr and their weak phases to be equal. Thus, 
we obtain the branching ratio in MSSM as 

Bt{B- ^ (/)n-)\MSSM < 1.2xl0-^ 

Br (S°^07r°) I M55M < 0.5 xlO-^ (22) 

Although the predicted branching ratios in MSSM are enhanced by one order from their SM 
values but they are still well below the present upper limits. 



4 Contribution from the VLDQ Model 

Now we consider the model with an additional vector like down quark [16]. It is a model 
with an extended quark sector. In addition to the three standard generation of quarks, there 
is an SU{2)l singlet of charge —1/3. The mixing of these singlet quarks with the three SM 
down type quarks provides a framework to study the deviation from unitarity constraint of 
3x3 CKM matrix. The important feature of this model is that it allows Z-mediated FCNC 
at the tree level. Thus, the presence of an additional singlet down quark imphes a 4 x 4 
matrix Via {i = 'W,c, t, 4, a = d,s,b,b'), diagonalizing the down quark mass matrix. For 
our purpose, the relevant information for the low energy physics is encoded in the extended 
mixing matrix. The charged currents are unchanged except that the Vckm is now the 3x4 
upper submatrix of V. However, the distinctive feature of this model is that FCNC enters 
the neutral current Lagrangian of the left handed downquarks : 



2 cos 9w 
with 



ULil^ULi - dLaUal3l^dL/3 - 2 siu^ OwJem ' (23) 



Uaf3 = E VliViP = - V:,V,p , (24) 

i=u,c,t 
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where U is the neutral current mixing matrix for the down sector which is given above. As 
V is not unitary, U 7^ 1. In particular its non-diagonal elements do not vanish : 

Uap^-V:^V^^^Q fora^/3. (25) 

Since the various Uap are non vanishing they would signal new physics and the presence of 
FCNC at the tree level, this can substantially enhance the branching ratio of S — > 07r. The 
observed discrepancy between StpKs S.^Ks can be explained in this model [17]. The new 
element Udb which is relevant to our study is given as 

Uab = KdKb + V:^V,, + V:,Vt, . (26) 

Thus the decay modes B ^ (f)7r receive the new contributions from the Z-mediated FCNC 
transitions and the new additional operator is given as 

(27) 

where Cy and C% are the vector and axial vector Zss couplings. Using the identity (Cf — 
Ca75) = [(Cf + C%){1 - 75) + (Cf - C%){1 + 75)]/2, the effective Hamiltonian for B ^ (pn 
transition is given as 



Kff - 



{C^ + CX){db)v-A{ss)v-A + {C^-CX){db)v-A{ss)v+A , (28) 



where the subscripts in the currents denote the usual {V — A) and {V + A) currents. Now 
evaluating the matrix elements of the operators the transition amplitude for the process in 
VLDQ model is given as 

AVLDQ^B- 4>^-) = ^U,,X{e-ps) ((Q + CX){1 + ^F^) + (Q - C^)) , (29) 

where we have also included the leading order nonfacorizable contributions. Now using the 
value for Cf and C% as 

C^ = ~ + lsm'ew , CX = -l, (30) 

sin^ ^VK=0.23, alongwith \Udb\ < 1-2 x lO'^ [18], we obtain the ratio of VLDQ and SM 
amplitudes as 

A^^^Q(E- ^ 07r-) 

The branching ratio in VLDQ model obtained using Eq. (16) as 

By{B' ^ cPn-)\vLDQ < 4.6x10-^ 

Bt{B' ^ cbn'>)\yLDQ < 2.1xl0-^ (32) 

which are two order above the standard model prediction and the branching ratio for B^ 
4>n^ is in agreement with the present data (1). 
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5 Contribution from the TC2 Model 



Now we calculate the branching ratio for B (pTC process in the framework of the topcolor 
assisted technicolor model (TC2) [19, 20]. It is well known that technicolor is one of the 
important candidates for the electroweak symmetry breaking and the extended technicolor 
model was proposed to generate the ordinary fermion masses. In order to generate large top 
quark mass the topcolor model has been constructed recently. Apart from some difference 
in group structure and/or particle contents, all TC2 models have similar common features. 
Following the TC2 model of Hill [19], the charmless decays B — > PP, PV are studied in Ref 
[21], using the generalized factorization approach. In this paper we will also use the same 
procedure to evaluate the new physics contribution to the B ^ (f)7i mode, but we will use 
QCD factorization approach to evaluate the hadronic matrix elements. 

In TC2 model, there exist top-pions (tt^ and 7f°), charged and neutral 6-pions {H^, 
and A'^ ) and technipions (vrf and tt^). The coupling of top-pions to t- and 6-quarks can be 
written as 

™* r ™* n 

(33) 



ittir'^ + +itRbLTT^ H -tihRTX^ + h.c 



ml 



where ml = {l—e)mt and ml ~1 GeV denote the masses of top and bottom quarks generated 
by the topcolor interactions and P^- is the top-pion decay constant. At low energy, potentially 
large FCNC arise when the quark fields are rotated from their weak eigenbasis to their mass 
eigenbasis, realized by the matrices Ul^r for the up-type quarks and D^ r for the down type 
quarks. Thus for example, making the replacement 

bR ^ D'^dR + D'isR + D'j^bR, (34) 

the FCNC interactions will be induced. In the TC2 model the corresponding flavor changing 
effective Yukawa couplings are 



{D'ltRSL + D^'^tRdL) + iH^ {D^'itLSR + D^ihdi^ + h.c. 



(35) 



The constraints on the parameters in the TC2 model are discussed in detail in Ref [21], 
obtained from various experimental data. For the mixing matrices the " square root ansatz" 
is considered, i.e. D'^^ = V^d/2 and D^f = l^ts/2. The other parameters are given as 

m^i = 100 GeV, m^g = 200 GeV, = 50 GeV, = 120 GeV e = 0.05 , (36) 

where P,r is the technipion decay constant. 

In this model, the decay process b — > dss is induced by the exchange of the charged top 
pions -n"^ and technipions nf and ng through the strong and electroweak penguin diagrams. 
Combining the new physics contributions with their SM counterparts, the effective Wilson 
coefficients are evaluated. 
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Thus the new strong and electroweak penguin diagrams can be obtained from the corre- 
sponding penguin diagrams in the SM by replacing the internal lines with the charged 
top-pions and technipions. The analytic expressions for these diagrams are calculated in the 
dimensional regularization with MS renormalization scheme in Ref. [21]. The new contribu- 
tions to the Wilson coefficients arising from these diagrams are given as 



r<TC2 
^0 



E' 



TC2 


<TC2 


TC2 



1 



1 

V2Gf 
1 

2V2Gf 



4F| 
1 

4^ 



3F2 



loiVt) + ^ ( loizt) + 8/o(et) + QNoi^t] 



SF3 



(37) 



where yt = nif' /m\ with mj' = (1 — e)mt, Zt = {trrit)'^ / m^^ and = {^''^h)^ /'>TT-'ig- The loop 
functions To (x), Fq{x), Io{x), Kq{x), Lo{x), No{x) are found from Ref [21]. Using the top 
pion and technipon masses to be 200 GeV, the values of the Co, Dq, - ■ ■ functions at the Mw 
mass scale are obtained as [21] 



{Co, Do, ,Eo, E'^y\=M^^{1.27, 0.27, 0.66,-1.58} 



(38) 



Now combining these values with the corresponding SM values at the H^-boson mass 
scale ( {Co, Dq, ,Eo, ^o}^^^U=Mv^ = {0.81,-0.48, 0.27, 0.19}) , and running the resulting 
contributions to the b quark mass scale (// = 2.5 GeV) the effective Wilson coefficients are 
obtained as [21] 

C3 = 0.0195, C4 = -0.0441, C5 = 0.0111 , Ce = -0.0535 , 

C7 = 0.0026 , Cg = 0.0018 , Cg = -0.0175 , Cio = 0.0049 . (39) 

Now substituting these values in (5) the branching ratios in the TC2 model are found to be 

Bt{B- 07r-)|TC2 = 1-2 X 10-^ , 

Br(S° ^ 07r°)|TC2 = 0.5 xlO-^ (40) 

Although the branching ratios are one order higher than the corresponding SM value but 
they are well below the present experimental limits. 



6 Conclusion 



In this paper, we have analyzed the decay mode B ^ (p-n both in the standard model and 
some beyond standard model scenarios. This decay mode proceeds through the quark level 
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FCNC transition b — > dss, receiving contributions only from one-loop penguin diagrams. 

However, because of OZI suppression, the QCD penguins play only a minor role in this case 
and the dominant contributions coming from the electroweak penguins. Therefore, in the 
standard model these decays are highly suppressed, which makes them a very sensitive probe 
for new physics. 

Using QCD factorization approach, we found the branching ratios in the SM for these 
modes as ~ O{10~^), which are quite below the present experimental upper limits O{10~^). 
We have also calculated the branching ratios in the MSSM with mass insertion approxima- 
tion, in the VLDQ model and in the topcolor assisted technicolor model. The branching 
ratios obtained in the MSSM and TC2 models are of the order of C(10~®), whereas they are 
found to be 0(10-^) for VLDQ model. 

Recently, BaBar Collaboration [4] has reported the first measurement of the branching 
ratio for — >• ^tt^ process as (0.21q;3 ±0.1) x 10~^. Although the error bars are quite 
large, but this preliminary experimental value is almost 2 orders larger than the standard 
model prediction. If in future the data will remain in this order, i.e., O{10~^), then it will 
give a clear signal of new physics effect present in the penguin induced process B (pir. As 
shown, only the VLDQ model can predict such a large branching ratio. Therefore, the future 
experimental data on B ^ (pn will serve as a very good hunting ground for the existence 
of new physics beyond the SM and also support/rule out some of the existing new physics 
models. 
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